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Abstract 

We develop a hydrodynamic approach to non-equilibrium conformal field theory. We study non¬ 
equilibrium steady states in the context of one-dimensional conformal field theory perturbed 
by the TT irrelevant operator. By direct quantum computation, we show, to first order in 
the coupling, that a relativistic hydrodynamic emerges, which is a simple modification of one¬ 
dimensional conformal fluids. We show that it describes the steady state and its approach, and 
we provide the main characteristics of the steady state, which lies between two shock waves. 
The velocities of these shocks are modified by the perturbation and equal the sound velocities of 
the asymptotic baths. Pushing further this approach, we are led to conjecture that the approach 
to the steady state is generically controlled by the power law t -1 / 2 , and that the widths of the 
shocks increase with time according to f 1 / 3 . 


1 Introduction 

One of the most powerful ideas in studying the dynamics of quantum field theory is that emerging 
from a hydrodynamic description of local averages |1]. Hydrodynamics allows to encode in 
a simple way non-equilibrium states, including states with constant flows and approaches to 
steady states, by concentrating only on quantities of physical relevance and without having to 
deal explicitly with an infinite number of degrees of freedom. In the quantum context, the 
passage from strongly interacting many-body quantum dynamics to classical hydrodynamics 
involve subtle effects that are sometimes hard to control [2]. It is therefore important to have 
simple models and setups where this passage can be studied with more precision. 

The main purpose of this paper is to implement the hydrodynamic approach in the context 
of perturbed non-equilibrium one-dimensional conformal field theory (Neq-CFT), and show that 
it emerges from the quantum description. 

Recall that Neq-CFT has been introduced to study out-of-equilibrium quantum phenom¬ 
ena using the extensive toolbox of conformal field theory (CFT) combined with an S-matrix 
approach [3]. It aims at describing non-perturbatively, and beyond the linear response theory, 
the low-energy sector of gapless one-dimensional quantum systems driven far from equilibrium. 
The driving is obtained by unitary evolution of independently thermalized subsystems, a setup 
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considered in other contexts before H El E]. Exact results in one-dimensional Neq-CFT [3] 
have been verified numerically in the Heisenberg chain |7]; other exact results along this line 
of thoughts have been obtained for free-fermionic quantum chains la ei mu hd na na and for 
higher-dimensional free models mug, and conjectures proposed for certain integrable systems 
mini- In higher-dimensional non-equilibrium CFT, a hydrodynamic approach was used within 
this setup in [I8][T9], leading to conjectured analytic formulas extending some of the results of 
[3] and numerically verified within the Gauge/Gravity duality [20]. As explained in [19], the lack 
of integrability plays an important role for the hydrodynamic argument, and, combined with 
other ideas, the hydrodynamic data give the exact quantum density matrix of the steady state. 

The CFT description of gapless systems is never exact but only asymptotically exact at 
low energy m- The asymptotic approach to the low energy effective dynamics is controlled 
by the so-called irrelevant operators - those which, by definition, do not influence the low- 
energy effective dynamics at equilibrium. Perturbations of CFTs by irrelevant operators are 
then important, especially far from equilibrium, in order to understand low-energy numerical 
[7] and experimental [22] results for non-equilibrium steady states and their universality [23 . 
[23], driven systems [25j and quantum quenches [22]. In this paper we study Neq-CFT with a 
particular irrelevant perturbation related to effects of band curvature; as this breaks integrability, 
hydrodynamic ideas are expected to be fruitful in this situation. 

2 A hydrodynamic approach 

Before embarking into the hydrodynamic approach, let us recall some facts about Neq-CFTs 
[3]. In the simplest case, one generates an out-of-equilibrium state in a CFT by first preparing 
two copies of this CFT on a semi-infinite half line each at different temperatures 7] = Pf 1 and 
T r = Pp 1 , and then gluing them at some initial time t = 0 through a contact point, say the 
origin x = 0. It was shown that a non-equilibrium steady state emerges on any interval [-L, L] 
for all times t > L/vp where vp is the Fermi velocity; beyond the interval, the system remains 
in its thermal state. That is, two shock waves are associated to this energy flow, which are sharp 
within the CFT description and which propagate at the Fermi velocity, along the two branches of 
the light-cone. The interior region of the light-cone is in a non-equilibrium state carrying a mean 
energy current Jp = ^ (T 2 — T r 2 ) with c the central of charge of the CFT, whereas the regions 
at the left and at the right of the light-cone remain at thermal equilibrium at temperature 7] 
and T r respectively. Since there are no fixed reservoirs in this setup, in general local energy 
gradients decrease as time evolves, so that the emerging steady state must be controlled by 
ballistic transport. A similar framework for generating non-equilibrium steady states was first 
studied in the context of electronic transport mm and in harmonic chains 0 and originally 
named the partitioning approach. 

Within the CFT approximation, the Hamiltonian is Hcft = f dx (T(x) + T(x)) where T(x) 
and T(x) are the two chiral components of the stress tensor [27]. They have scaling dimension 
2, so that the Hamiltonian has the correct dimension 10. The most natural irrelevant operators 
are bilinear in the stress tensor components, that is TT, and T 2 or T 2 . For simplicity, we 
shall restrict ourselves to the TT perturbation, associated to curvature effects in the dispersion 
relation of the underlying gapless microscopic model. The perturbed Hamiltonian we shall 

1 By convention, we normalize the Fermi velocity to vf = 1. We set ks = 1 and h = 1. 
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consider is thus 


H = dx (T(s) + T(x)) + g / d xT(x)T(x) 


' 


(1) 


At equilibrium this perturbation is irrelevant because the operator TT has dimension 4 and 
hence g scales like a? with a the short distance (UV) scale of the microscopic model. 

In the following sub-sections we use a hydrodynamic approach to study the effects of the TT 
perturbation on the non equilibrium steady state, on the approach to the steady regime, and on 
the structure of the associated shock waves, and we show by a direct perturbative computation 
that these hydrodynamic considerations are correct. 


2.1 Generalities: hydrodynamic description near criticality 


The hydrodynamic approach is rooted in two principles [2S]: the first consists in assuming that 
there is some kind of local thermalization and the second that the equations of motion follow 
from conservation laws. 

Local thermalization assumes that at every point x and time t, the state may be described, 
locally in a neighborhood TX X) t,. by an effective density matrix pi oca i of the form of an equilibrium 
density matrix, but involving all available local conserved charges of the theory. Near criticality 
in non-integrable models, omitting internal symmetries as we are only considering energy trans¬ 
port, the only charges that may describe transport are the hamiltonian H and the momentum 
P, associated to invariance under time and space translations. In non-integrable systems there 
are no other local conserved charges. With h and p the hamiltonian and momentum densities 
respectively, piocal is of the form 


piocal oc exp 


/ 3h(x,t ) / dx' h{x') + /3 p (x, t) / dx'p{x')\ 

J%*.t -Arw J 


( 2 ) 


with x, f-dependent parameters (local potentials) (3h and f3 p . This implies that the expectations 
of local observables can be computed as in a generalized thermal state of the form e~P hH +Pp p 
(generalized by the presence of the momentum operator P), but with local effective inverse 
temperature /3h and momentum potential f3 p that depend on space and time. 

In quantum chain models, there is of course no continuous space translation, hence no con¬ 
served momentum operator. Continuous space translation only emerges near criticality, where 
low-energy excitations are supported on such large scales that the underlying lattice structure 
is not important. Without a conserved momentum operator, the density matrix is locally ther¬ 
mal, hence it cannot describe energy transport^. It is the presence of a conserved momentum 
near criticality that allows, within this hydrodynamic description, for a possibly nonzero energy 
current to develop at large times in the partitioning approach. This agrees with the intuition ac¬ 
cording to which only ballistic currents may develop in this approach, while the lattice structure 
introduces diffusion away from criticality. 

The choice of the local observables that should be described by local thermalization is im¬ 
portant. For instance, in a thermal state, averages of derivatives of local observables are zero by 
translation invariance in space and time. On the other hand, in a locally thermalized state the 
derivative of the average of an observable is generically nonzero. Hence in a locally thermalized 


2 With integrability, this problem can be circumvented by the presence of appropriate non-trivial conserved 
charges. 


3 



state, the choice of putting the derivative inside the average or outside matters, and this is 
guided by the physical scale on which the derivative is to be taken: either at the microscopic 
level (inside a hydrodynamic cell, where quantities are effectively constant), or at the hydrody¬ 
namic level (there are variations from one hydrodynamic cell to the next). In the hydrodynamic 
description of QFT one assumes that conserved densities and currents are locally thermalized, 
and that hydrodynamic equations are determined by the conservation laws. 

In the present situation we have two local conservation laws, those of time and space trans¬ 
lation invariance: 


d t h + d x j = 0, d t p + d x k = 0. (3) 

Within the (pure) hydrodynamic approximation, expectations of the densities h and p and the 
currents j and k are evaluated by local thermalization: they are first evaluated in generalized 
thermal states, then the local potentials are made dependent on space and time. Hence they 
become functions of the local potentials f3h(x,t) and f3 p (x,t). Let us denote these expectations 
by h, p for the densities and j and k for the currents. From eq. (J3|) we get two equations 

d t h + d x j = 0, d t p + d x k = 0. (4) 

These are non-linear differential equations for fih and /3 P , which are the hydrodynamic equations. 
Notice that we have two equations for two functions with given initial condition, so that in 
principle the problem is reduced to solving non-linear PDEs. 

In thermal equilibrium, the energy current and momentum density are zero by time-reversal 
and space-reversal symmetry. Since there is only one parameter (the temperature), there is an 
equation of state k = F( h) relating the pressure k to the energy density h. In the generalized 
state, the equations of state may be generalized to k = F(h,p), and j = G(h,p). With the 
generalized equations of state, the hydrodynamic problem is recast into equations for local 
energy and momentum densities, instead of equations for /3h and /3 P 

Other parametrizations of the hydrodynamic problem are possible, and one often uses “fluid 
velocities”. In some situations, it is possible to use a parametrization that explicitly separates 
“dynamical” degrees of freedom from “thermal” ones. Assume that there is an additional dy¬ 
namical symmetry of boost: a space-time transformation operator B that transforms energy into 
momentum, [B,H] = P, and that preserves the space spanned by H and P. Two immediate 
examples are the relativistic and the galilean boosts, [B,P] = H (with Casimir H 2 — P 2 ) and 
[B,P\ = ml (with Casimir H — P 2 /(2m)), respectively. The former applies to critical points 
with dynamical exponent z = 1, the latter with z = 2 (see for instance [29] for a hydrodynamic 
study of non-equilibrium states in this context). Using the boost operator, we can always write 
(3hH — (3 p P = /3 res t e~ eB He dB . The choice of B is not unique, and for every space-time point 
x,t we have a boost B x j that keeps x,t invariant. The local density matrix can be written as a 
local boost of a local thermal state, 

Piocal OC exp ( - /3 rest (x, t) [ dx' e ~ e ( x ’Q Bx ’ t h(x')e 0 ( x ’ 1 ) Bx A . (5) 

V J%c,t J 

The local boost parameter 9(x,t) is the dynamical part, and the local rest-frame inverse tem¬ 
perature /3 res t(x,t) is the thermal part. This hydrodynamic parametrization was used in [19] 
in order to give a characterization of higher-dimensional non-equilibrium CFT and describe the 
steady state as a boosted thermal state. 
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It turns out that the requirement that the momentum density be equal to the energy current 
up to, possibly, an even derivative term, p = j + p d'^j , is essentially enough to imply relativistic 
invariance. Further, this implies that, in a thermal state, 

= h + k , (6) 

and that, given the thermal equation of state k = F( h), all generalized thermal averages are 
explicitly expressed as functions of 6 and h, see eq. (!54p in Appendix [A] These facts, which we 
show within the general context of local many-body systems in Appendix [Al can be useful to 
simplify the description of the hydrodynamic problem. 

Finally, the hydrodynamic approximation discussed above is the “pure” hydrodynamics. It 
can be made more accurate by adding derivative corrections, including viscosity terms. These 
are corrections, in the expressions for the densities and currents, involving higher derivatives of 
the basic quantities (for instance corrections to the equations of state involving derivatives of 
energy and momentum densities), and are associated with shrinking the hydrodynamic cell in 
order to take into account more and more of the microscopic variations. 

We will make all these concepts more precise in the example of the perturbed CFT studied 
in the present paper. 

2.2 The exact first-order hydrodynamic description of TT perturbed CFT 

In the present setting of perturbed conformal field theory, the hamiltonian density could be taken 
as T(x) + T(x) + gT(x)T(x), and the momentum density p(x) = T(x) —T(x). Of course, these 
are defined up to total derivatives of local densities, and for our purposes it will be convenient to 
add the term —gc/{24'K)d 2 (T{x) + T(x)) to the above energy density. The equations of motion 
are then determined by the Hamiltonian (jT]) , via dtO = i[H,0] for any operator O. The basic 
commutation relations for the stress-energy tensor are 

-i[T(x),T(y)\ = -(T(x) + T(y )) 5'(x - y) + (x - y), 

with c the central charge of the CFT m, and similarly for the T commutation relation except 
for the change of i —>■ —i. The T and T components commute. Hence, we have 

h{x) = T{x) + T(x)~-^(T"{x) + T"{x)) + gT(x)T{x) 

z47T 

j(x) = T(x)-T(x) + 0(g 2 ) (7) 

for the energy density and current, and 

p(x) = T(x) — T(x) 

k(x) = T(x)+f(x)-^(T"(x)+T"(x))+3gT(x)T(x)+0(g 2 ) (8) 

for the momentum density and current. Here and below primes denote space derivatives. 

Let us briefly analyze the case g = 0. In this case, the local density matrix is 

Piocal oc exp ( - [ dx'T(x') - (3(x,t) [ d x'T(x')\ (9) 

v -Fn x , t J^ x ,t J 
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with x, t-dependent parameters (3 = f3h — /3 P and (3 = /3h + /3 P , and we have j = p = T — T 
and k = h = T + T. Hence the averages of T and T can be used as hydrodynamic variables, 
and the expectations of any local observables, say products or powers of T and T at some point 
x and time t, can be computed, within the hydrodynamic approximation, as in equilibrium 
CFT but with local effective “temperatures” (3(x,t) and /3(x,t). In particular, there is chiral 
factorization: right-movers and left-movers are independently thermalized. At g = 0 the CFT 
enjoys relativistic invariance, with boost operator B = —i f dxxh(x) satisfying [B,H] = P and 
[B, P ] = P[, and the above can be interpreted as a boosted density matrix with rest-frame 
inverse temperature /3 res t = \/]3]3 and boost parameter given by tanh# = (]3 — f3)/(]3 + (3). The 
thermal average at temperature r is [301 El] 


(T) = <T> = -r 2 


( 10 ) 


and with the initial condition where left and right halves are independently thermalized at 
temperatures Tj and T r respectively, the hydrodynamic problem gives (3(x,t ) = T^ 1 @(x + t ) 
and f3(x,t ) = T r -1 0 (—x + t). As was shown in [3, IT)) . this hydrodynamic description is in fact 
exact, and does not necessitate any derivative corrections: it is a direct consequence of chiral 
factorization at the level of the quantum equations of motion. This solution represents sharp 
shocks emanating from the contact point and propagating at the speed of light in opposite 
directions, between which the steady states lies. 

Let us now consider the first order correction in g. We first use the facts that, in a generalized 
thermal state, there is translation invariance, (T") = (T") = 0, and that to leading order in 
g there is chiral factorization g(TT) = g(T)(T) + 0(g 2 ). We can then parametrize averages 
using the two quantities w(x,t ) = ( T(x,t )) and w(x,t) = (T(x,t )) evaluated within the pure 
hydrodynamic approximation using pi OC al) and we get 


h = w + w + g ww + 0(g 2 ), j =w — w + 0(g 2 ), 
p = w — w, k = w + w + 3g ww + 0(g 2 ). 


( 11 ) 


Recall that the equations of motion are dth+d x j = 0 and Ojp+clxk = 0. These are two equations 
for the two unknowns w and w, chosen as a parametrization of the two unknown local potentials 
f3 h and f3 p . 

This description is not exact at order g: there are derivative corrections. Fortunately, in the 
present case it is possible to obtain the full, exact derivative corrections at order g: they are 
a direct consequence of the quantum equations of motion. It turns out that these derivative 
corrections are “trivial”: they can be entirely absorbed into appropriate definitions of basic 
hydrodynamic quantities, in such a way that the hydrodynamic equations remain unchanged. 
Writing the conservation laws for quantum averages using (0 and ([8]) without assuming local 
thermalization, using only chiral factorization at g = 0, the definitions 


w(x,t) = (T{x,t )) - — «T"(x,t)) + (T"(x,t))) 
w{x,t) = ( T{x,t)) - (( T"(x,t )) + (' T"(x,t ))) 


( 12 ) 


give rise to the pure-hydrodynamic equations for the quantities (HD. That is, in terms of these 
exact quantum averages, the conservation equations for the densities and currents (HID are a 
consequence of the quantum equations of motion, hence are exact for all x and t to order g. All 


6 


derivative corrections have been absorbed into the definitions ( 1121 ) and the choice of the energy 
density in (|7|). 

The initial conditions are those reflecting that initially the left and right parts are inde¬ 
pendently thermalized at temperatures T) and T r respectively. At thermal equilibrium with 
temperature r, the average momentum density is zero (and averages are space-time indepen¬ 
dent), hence w = w =: w t h(r). By dimensional analysis 

w th{ T ) = j^T 2 f{gT 2 ) (13) 

for some function / with /(0) = 1. We expect / to have an expansion in positive integer powers 
of g, and the normalization of the zeroth order m gives 

f(u) = 1 + au + 0(u 2 ). (14) 


Hence, at initial time, 


w = w = —I} 2 (l + agT 2 + 0(g 2 )) on the left (x < 0) 
w = w = ^T 2 (l + agT 2 + 0(g 2 )) on the right (x > 0) (15) 


with an abrupt crossover in the neighborhood of the contact point (that is, close to the origin). 
Below we compute the constant a by using relativistic invariance, giving 

C7T 
6 ~ 

This fully fixes the hydrodynamic problem (with all derivative corrections) at order g. 



Remark 2.1 Given any solution w(x, t ) and w(x, t ), the quantum averages {T(x, t)} and {T(x, t)), 
obtained by solving m, are naively not uniquely fixed: 

(T(x, t)) = w(x, t ) + (w"(x, t ) + w"(x, t)) + uj(x, t ) 

4<37T 

(T(x, t)) = w(x, t ) + (w"(x, t ) + w"(x, t)) + u>(x, t ) 

4<37T 

where ui(x,t) is solution of the equation ui — = 0(g 2 ) whose kernel is of the form 

A + e a: /v / 9 c /( 247r ) + A_ e~ x !Vd c /( 247r ). However, this kernel is non-perturbative in g, hence beyond 
the present order-g calculation. 


2.3 Equation of state and relativistic structure at order g 

The thermal equation of state relating the pressure k to the energy density h can be obtained 
from (HID by using the fact that w = w in a thermal state: 

k = F(h) = h+|h 2 + 0( 5 2 ). (17) 

The generalized equations of state can be obtained similarly, 

k = h + |(h 2 - p 2 ) + 0(g 2 ), j = p + 0(g 2 ). (18) 
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Associated to the equilibrium equation of state m is the sound velocity (at temperature r) 


v s = V F'{ h) = 1 + |h + 0(g 2 ) = 1 + ^r 2 + 0(g 2 ) (19) 


Remark from (HID that to order g, we have j = p. This equality implies propagation of small 
disturbances of energy, momentum and pressure densities near their equilibrium values occurs 
through sound-like waves with this sound velocity (see for instance the arguments presented in 
[321 [33]). We will see below that the sound velocities also play a role far from equilibrium at 
order g. 

Further, observe from (ED and ([HD that the equality j = p holds at the operator level at 
order g. This implies that relativistic invariance emerges (see Appendix [A] for the general 
theory). Indeed one can verify explicitly that the boost operator B = —i f dxxh(x), with h{x) 
the perturbed energy density ED) satisfies [B,H] = P + 0(g 2 ) and [B,P] = H. Therefore, 
the averages in the generalized thermal state are expressed in terms of relativistically boosted 
observables in a thermal state. Using the result (1541) along with (1171) . we may express them as 


h = cosh 29 h rest + ^ (cosh 29 - 1) h 2 est + O (g 2 ) 
j = P = sinh 29 (h rest + | h 2 est ) + 0{g 2 ) 

k = cosh20h rest + |(cosh20 + l)h 2 est + 0(g 2 ). (20) 


where h rcgt is the rest-frame thermal energy density h rest = h t h(T res t), deduced from (flTT) . (fl3l) 
and CED, given by h th (r) = ^r 2 (l + gr 2 (a + + 0(g 2 )). The relation Q can be combined 

with the thermal equation of state EZD and the form (1211) of the thermal energy density. A 
simple calculation shows that this fixes a to (fl6D . and thus 


hth(r) 


C7T 9 

-7 

6 



( 21 ) 


Remark 2.2 Formally extending the above thermodynamic equations to large values of |(?| one 
finds that, in some respects, g < 0 appears to be more natural. For instance, the energy density 
EID has an unphysical maximum, as a function of temperature, if g > 0. Further, if g > 
0, the sound velocity (ED is greater than the Lorentz speed of light (which is set to unity) 
associated to the (order-g) relativistic invariance. However, these facts do not lead to perturbative 
inconsistencies, and none of the perturbative calculations below are affected by the sign of g. 
The temperature at which the energy density is maximal is nonperturbative (of order 1 /g); and 
a Lorentz transformation from a point in the space-time region lying between the light cone and 
the sound cone to a point at a time t = 0 would involve a very large boost parameter (non¬ 
perturbative in g), thus precluding, in perturbative considerations, the conclusion of nonzero 
equal-time correlations (breaking of causality). In fact, the Lieb-Robinson velocity vlr can be 
shown perturbatively to be greater than 1 + Ag for any constant A > 0, and thus greater than 
the sound velocity no matter the sign of g. Below we keep g arbitrary, and all conclusions stay 
valid perturbatively in g independently of its sign. 



2.4 Light cone effects, shocks and steady states 

In this section, we solve the hydrodynamic problem (J4|) associated to the TT perturbation, using 
the parametrization in terms of the functions w and w in Eqs. (ED- We decipher how the TT 
perturbation modifies the shock propagation, the structure of the shocks as well as the approach 
to the steady states and the steady state itself. We shall work to first order in g only. 

2.4.1 Shock wave assumption 

The simplest analysis of the hydrodynamic problem is that which assumes that the picture of 
shocks emanating from the contact point, which holds exactly at g = 0, is accurate even with 
g / 0 at large scales. This picture was used in HSUS] in order to generalize results to higher¬ 
dimensional CFT. More precisely, we make the assumption that, with exponential accuracy, 
averages of local observables as functions of the space-time point (x, t ) take steady and uniform 
values whenever |x + vit\ 2> £/(f) and \x — v r t\ £ r (t), for some £i, r (t) bounded from above by 
a power law oc t a for all t with some exponent a < 1. This represents shocks at speeds vi > 0 
and v r > 0, propagating towards the left and the right respectively, with time-dependent widths 
£ i(t ) and £ r (t), respectively, which grow sublinearly at large times. By the initial conditions, 
two the left (right) of both shocks, averages are evaluated in the initial left (right) reservoir, and 
the steady state lies between both shocks, described by a generalized thermal state. 

It turns out that this shock assumption, along with the knowledge of the averages in the 
generalized thermal states, is sufficient to fix both the shock velocities and the steady-state 
parameters. It does not, however, give any information on the approach to the steady state and 
on the structure of the shocks themselves. 

Consider the conservation laws ((3|) integrated along a rectangular contour crossing a shock 
and whose diagonal goes from (xi,fi) to (x 2 ,t 2 )- That is, 

rx 2 rt 2 

/ dx (h(x,t 2 ) - h(x,ti)) + / dt (j(x 2 ,t) - j(x ly t)) = 0, 

Jx 1 J t\ 

with (xi, ti) and (x 2 , £ 2 ) both lying inside the shock, and similarly for p and k. These equations 
can be written for the left and the right shocks. We can choose the rectangle such that |x 2 — xi| = 
Vi )T \t 2 — ii| with |x 2 — xi| much larger than the width of the shock. Using the assumptions that 
the densities and the currents are asymptotically uniform and steady away from the shocks, that 
the shocks are of extent that grows sublinearly in time, and that the current vanishes outside 
the light-cone, this yields two sets of two equations, one for the left shock the other for the right 
shock, 

vi(hi - h s ) = j s , vip s =ki-k s , v r (h s - h r ) = j s , v r p s = k s - k r , (22) 

where indices l, r, s indicate the averages in the left reservoir, right reservoir and in the steady- 
state region, respectively. Adding them allows to eliminate h s and k s , keeping only the densities 
h i <r and the pressures kevaluated in the asymptotic region at equilibrium, 

(vi + v r ) p s = k 1 - k r , (u z _1 + u" 1 ) j s = hi - h r . (23) 

Combining them, we find the two equations 

(hi - h s )(k 1 - k s ) = (h r - h s )(k r - k s ) = j s p s . (24) 
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These are two equations for two steady-state unknowns. The steady-state unknowns can be 
taken as w s and w s using the parametrization (1111) , the boost and rest-frame temperature 6 and 
T rest using the relativistic parametrization (12U1) . or simply the steady-state energy density and 
current h s and j s using the generalized equations of state (1181) . A calculation up to order g 
gives, within the latter parametrization, 

h s = + 2 hr + |(h i -h r f + 0{g 2 ) 

js = 2 r + - K) + 0{g 2 ). (25) 

The thermal averages h; r can be evaluated in terms of the temperatures T) r using (1211) . and in 
particular we obtain 


Js = 


C7T n 

= —T? 


12 


(l - 9 -^t 2 ) 

C7 T 9 

- T ? 

f 1 - 9 -^t 2 ) 

V 12 1 J 

12 r 

V 12 r J 


+ 0{g 2 


The shock velocities can now be evaluated through 


VLr = 


h l : r h s 


= 1 + 


9 ^^ 2 
12 


Ttr + 0(g 


(26) 


(27) 


We make several observations. 


• First, to this order in g, the steady-state current is still a difference of a function of 
the left-reservoir temperature minus the same function of the right-reservoir temperature, 
j s = J(Ti) — J(T r ). In the unperturbed case this had important echoes on the full counting 
statistics, and in general this implies that the non-equilibrium current can be obtained 
purely from the linear-response conductivity G(r) = dj s /dTi\ T=T _ t = dJ(r)/dr as j s = 

It',. dr G ( t ) 0 - 

• Second, the shock velocities vi and v r are exactly the sound velocities (flfll) of the left and 
right reservoir, respectively. That is, it is the linear sound velocities of the reservoirs that 
control the speed of the shocks describing the far-from-equilibrium steady state. 

• Third, parametrizing the reservoirs using w (with the thermal result (1131) . (1141) . (1161) 1 we 
note that the current takes the simple form 

js = viwi - v r w r + 0(g 2 ). (28) 

This can naturally be interpreted as the ballistic transport of energy by independent 
left/right movers, with chiral energy densities wi r and velocities vj r that are solely deter¬ 
mined by the reservoirs’ temperatures. One may understand this as a generalization of the 
case g = 0, where w and w are right- and left-moving quantities at the speed 1: to order 
g , w and w may still be seen as independent right- and left-moving energy densities, but 
their velocity is “dressed” into the sound velocity by the reservoirs’ state. This is also an 
extension of the picture that holds near equilibrium. With 7] ~ T r , there is linear wave 
propagation with the sound velocity v. The natural independent right- and left-moving 
energy densities are the combinations of space-time dependent small variations (near the 
equilibrium values) Se± = (5h ± u _1 (5j)/2, and the current is v(5e + — <5e_). 
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• Fourth, near equilibrium (i.e. with 7} ~ T r ), linear wave propagation implies that [33| 
js ~ (k; — k r )/(2u) where v is the equilibrium sound velocity. One can indeed check that 
this agrees with (1281) . using d(vw)/dr = (l/2v) dk/dr (where v, w, k are at equilibrium 
with temperature r). 

• Fifth, one can verify that the inequality is > (k; — k r )/(2 vlr) |33j is verified, where vlr 
is the Lieb-Robinson velocity (see Remark 12.21 concerning vlr). 

• Finally, we quote the steady-state rest-frame temperature and boost velocity from its 
relativistic parametrization: 

T rest = s/Wr (1 - - T r f + 0{g 2 )) , tanh 6 = (l - 

In particular, this means that in the steady state, the potential /3h 
hamiltonian in the density matrix as in Q takes the simple form 

a A + A- 

2.4.2 Solution to the hydrodynamic problem: appearance of shocks 

In the previous paragraph, we have made the assumption that two shocks emanate from the 
connection space-time point. In the present paragraph, we provide the solution to the full order-g 
hydrodynamic problem. The main observations from the calculation below are as follows: 

• The problem does not have a unique solution. The space of solution includes the two-shock 
solution above, but also contains solutions with additional “remnant” shocks at x = ±f. 

• Nevertheless, the steady state is unique: every solution gives rise to the same state in the 
central region (this state is in agreement with that found from the two-shock assumption). 

Let us discuss briefly the first point above. The non-uniqueness of the weak solutions to 
the Riemann problem (the step initial conditions) for pure hydrodynamic equations has been 
observed before m- In standard hydrodynamic problems, it is usually possible to lift this non¬ 
uniqueness by considering the entropy current or any other inequality condition [36]. With the 
requirement that the entropy cannot decrease at the shocks (a local implementation of the 2 nd 
law of thermodynamics), this may forbid certain shocks and lead to rarefaction waves (transition 
regions of extent growing linearly with t). Positivity of local entropy production occurs with 
viscosity terms, and weak solutions to pure hydrodynamics should be seen as emerging, at 
large scales, from viscous hydrodynamics; thus positive local entropy production at the shocks 
becomes a condition for selecting weak pure hydrodynamics solutions. In the present case, 
higher-derivative corrections are completely absorbed into redefinitions, hence at order g there 
is no positivity condition that can fix the solution. We hope to come back to this problem in 
the future by analyzing higher orders in g. 

Let us proceed with the calculation. The full order-g hydrodynamic problem is given, in 
terms of the parametrization w, w as defined in (1121) . by the conservation equations (J4]) with 
ED, along with the initial conditions (1151) . 


(30) 


^ Tl T r + 0(g 2 )). 

(29) 

associated with the 
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At g = 0, denoting w = wo, these equations simplify to d x wo = —dfWo and d x wo = dtwo- 
With the initial conditions, the solution is immediately given by 

= if ( T ' 20 ( t “ x ) + T r Q ( x ~ t)) i g -x) + T 2 Q(x + t)) . (31) 

This indeed represents two sharp shock waves both at velocity 1. 

Up to first order in g, we may solve the hydrodynamic problem by using this zeroth-order 
solution in order to simplify the bilinar terms proportional to ww in (HD. The main idea is that, 
at order g , we may replace gw and gw by gwo and gwo respectively. We first divide space-time 
(with positive times t > 0) into three regions: the left (L) x < t, the center (C) —t<x<t 
and the right (R) x > t. Within each of these regions, wq and wq are both constant, with 
wq = T) 2 in the center and left, wo = f§T 2 in the right; and wq = T f 2 in the center and 

right, w 0 = ffR/ 2 in the left. Hence we may write gdt(ww) ~ gwodtw + gwodtw with wo and wq 
piecewise constant, and similarly for gd x {ww). These transform the bilinear terms into linear 
terms with piecewise constant coefficients. The equations then become 

(d x + (1 + gw 0 )d t )w - (< 9 X - (1 + gw 0 )d t )w = 0 
(dt + (1 + 2>gw 0 )d x )w - (d t - (1 + 3gw 0 )d x )w = 0. 


Eqs. (1321) are linear equations which we can solve, within each region, by looking for solu¬ 
tions of the form e l (p x - E P t ) _ A straightforward calculation shows that these plane waves have 
dispersion relations with two branches: 


E _ [p{^ + gwo + 0(g 2 )) 

P { -pi 1 + gwo + 0{g 2 )). 


We observe that the dispersion relation, to order g, is still linear, and hence only changes the 
propagation velocities. This means that the solutions for w(x,t) and w(x,t ) have the form, up 
to 0(g 2 ) terms, 


w(x, t) = w + [x — v + t) + w {x + v t ), w{x, t) = w + {x — v + t) + w (x + v t ) (34) 

where v + = 1 + gw o and v~ = 1 + gw o- In the left (L), center (C) and right (R) regions we 
obtain the velocities: 


= v c = 1 + 


fl'CTT 2 


12 


T, 


l ’ 


4 


= v+ = l + 


gcir 2 
12 r: 


(35) 


and within these regions we denote the left- and right-moving waves by w± q r and c R . Note 
that the propagation velocities (1351) are the shock velocities vi and v r obtained in paragraph 12. 4. II 
from assuming the presence of two shocks. 

For definiteness we assume that g > 0, so that all velocities are greater than one; a similar 
calculation can be done for g < 0, leading to the same conclusions. 

The functions w + , w ~, w + and w~ in (1341) are related to each other. This can be obtained 
by putting (1341) into any one of the two equations (132[) . for instance the first one. Since the 
arguments of the functions are either x — v + t or x + v~t and the equation holds for all x and 
t (within a given region), this single equation leads to two relations, one for the right movers 
the other for the left movers. These two relations contain single-derivative terms with constant 
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coefficients, as in (132ft . Integrating, one obtains similar relations but without derivatives. The 
integration constants may be set to zero, which accounts for the “gauge symmetry” in (|34|) . 
w± i —y w ^ ± A and w± H > vu^ ± B for constants A and B. After these operations, we find that 
the relations are 

(1 ^ v + v ± )w ± = (1 ± v + v Zf )w ± + 0(g 2 ). (36) 

Since the velocities are of the form v± = 1 + O(g), we see that w + = O(g) and w~ = O(g). It 
is then convenient to work solely with w + and w~ instead, and using (|36ft and the form of 
we find that m becomes 

w = w + — gwQW~ + 0(g 2 ), w = w~ — gujQW + + 0(g 2 ). (37) 

The boundary conditions are conditions on the functions wr^l and wr^l- On the right, we 
have wr(x, 0) = wr(x, 0) = w r for all x > 0, and on the left we have wl(x, 0) = wl(x, 0) = wi for 
all x < 0, where we recall that are the thermal values Using the decomposition (1571) 

and the velocities (|35j) . as well as the leading part of the thermal values wi >r = (cn/12)T^ r +0(g), 
this gives 

w+(x) = iDx(x) =w r + gw 2 + 0(g 2 ) [x > 0) 

■ 0 0 ^oo J 

Wl(x) = u>l(x) = wi+gw l +0(g ) (x < 0 ). 

The full solution to the problem is obtained by supplementing the boundary conditions 
with the continuity conditions between the regions, at the interfaces x = ±t. The continuity 
conditions are obtained by considering the integrated version of the equations (1321) . which, we 
recall, are the full order-g equations and are equivalent to (0]). Consider first the jumps at the 
interface x = t between the center and right regions. For a shock of velocity 1 between the 
center and right regions across the interface x = t. this is (using j = p + 0{g 2 )) 

he - h R = j c ~ 3R = k c - kR. 


From (fill) , the first equation is equivalent to 

wr-wc = | (wc - wr)wr + 0(g 2 ), 

while the second is similar but with g replaced by 3 g. Combining, one possible solution is 
wr = O(g) ; but it turns out that this is inconsistent with the boundary condition. Hence the 
only consistent solution is 

wr-wc = 0(g), wr-w c = 0(g 2 ) (at x = t). (39) 

That is, to order g there is no jump in w at that interface, but there may be an 0(g) jump in 
w. A similar calculation at the interface x = — t gives 

w L -w c = 0(g ), w L -wc = 0(g 2 ) (at x = -t). (40) 

Let us now combine the continuity condition at the interface x = t with the boundary 
conditions. Using the decomposition (1371) and the velocities (135]) . the second relation of the 
continuity condition (1551) is 

Wq(( 1 + Vi)t) - gw r w£(( 1 - v r )t) = Wr((1 + v r )t) - gw r w^(( 1 - V r )t) + 0(g 2 ) 
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( 41 ) 


while the first is w^{{\ — v r )t) = w~^(( 1 — v r )t) + O(g), thus leading to 

Wc((l + vi)t ) = «£((1 + v r )t) + 0(g 2 ). 

With the initial conditions (13811 we then find 

u>q(x) = w r + gw 2 + 0(g 2 ) (x > 0). (42) 

A similar calculation combining the initial conditions with the continuity conditions (14011 at the 
interface x = — t gives 

Wc(x ) = wi+ gw 2 + 0(g 2 ) (x < 0). (43) 

These then imply, using the above assumption that the velocities are greater than 1, 

w\(x) = wi + 0(g) (x < 0), w2(x) = w r + 0{g) (x > 0). (44) 

Equations (|38jl fully fix wr to the right of the right shock, x > v r t, and wl to the left of the 
left shock, x < —vit. Further, since the velocities are greater than 1, equations (|42l) and (1431) 
fully fix wc in the center region, —t<x<t. The remaining parts are the areas between the 
regions’ interfaces and the shocks, t < x < v r t and —vit < x < —t. These are fixed by equations 
PUP , but only up to O(g) terms. Hence, the order g problem does have a unique solution. The 

space of solutions include the two-shock solution considered in paragraph 12.4.11 which is that 

obtained by setting w^(x) = wi + gwf + 0(g 2 ) (x < 0) and w~^{x) = w r + gw 2 + 0(g 2 ) (x > 0). 
But it also includes a continuum of solutions, which in general present four shocks, at x = v r t, 
x = —vit and also at x = ±t. 

Despite the non-uniqueness of the solution to order g, it is clear that the steady state is 
unique. Indeed, the steady state is completely determined by the center region, fixed by equa¬ 
tions (USD and (1431) and the decomposition (1371) . We find 

w c = wi- gwi(w r - wi), w c = w r + gw r (w r - wi), (45) 

and this leads to the following steady-state energy current and density: 

j s = w C -w c = (wi-w r )(l + g[wi + w r ) + 0(g 2 )), (46) 

h s = w c + w c + gwcw c = wi+w r + g(wf - wiw r + wj) + 0(g 2 ), (47) 

in agreement with the two-shock solution (1251) . 

2.4.3 Approach to the steady state and spreading of the shocks 

We now consider how the steady states is approached, and study the internal structure of the 
shocks. In order to do so, we extend beyond the first order in g the analysis by using properties 
we expect to be valid at higher orders. The following results are conjectural. 

Approach to the steady state. The steady state emerges at large time in the center region. 
In order to understand how the system approaches steadiness we have to look at the large time 
behavior of w and w in this center region, say at x ~ 0. Let us solve (1321) to all order in g. 
Continuing the solution in the central region all the way beyond the region, at t = 0 we find an 
initial condition that is a sharp profile (modified by the jumps through the shocks). Because of 
this sharp profile, it is better to consider the space derivatives w' and w ', whose initial conditions 
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are then of the form w'(x,t = 0) = w'(x,t = 0) oc 5(x). These have smooth Fourier transforms, 
and the solution can be written as 


w'(x « 0, t) ~ 



with Ep the two branches of the dispersion relation in the center region and U p and V p smooth 
functions of p. In the bulk of the steady state, we assume that we can solve recursively, order 
by order in g, the system of differential equations by bringing it to linear equations. Thus the 
above form should also be valid, but with modified dispersion relation. This is our working 
hypothesis. Assuming that U p and V p stay smooth (because we are dealing with the derivative 
of w), for x fixed the large time behavior of the integral representing w(x,t ) can be evaluated 
by a saddle point approximation. The saddle is at some p* and the saddle point integral is then 
of the form f dpe lconst ( p ~ p *> t . As a consequence, w'(x,t ) decreases as 1 /y/t. This behavior has 
been seen in the approach to the steady state of the non-equilibrium Ising model [T 2 j and in 
Ising quantum quenches [34 1. 

Estimating the width of the shocks. We consider a similar type of calculation in order to 
estimate the width of the shocks. For this purpose, let us assume that E± is of the form 
E± = ±p(a ± + 7 =l= p 2 -|-- • •), where the coefficients 7 ^ are of order g 2 . Although parity constrains 
the possible dispersion relations, we do not know how to justify within this picture the absence 
of a p 2 term in full generality, except that it is valid in the Ising case, but we will see that the 
result agrees with an independent argument based on quantum flows. 

We want to probe points at the “shoulders” of the shocks, so that we set xEv^t = et and we 
look at the large time behavior with e fixed but large enough. Again, we perform a saddle point 
approximation. A simple calculation gives the saddle point to be p* ~ y/e, and the saddle phase 


is e 


i conste 3 / 2 t 


Replacing et by x ± vf,t leads to the phase exp 


i const 


/ (x±vft) 3 


1 / 2 ' 


and thus 


we estimate the width of the shocks to be 5£ s hock — (gi) 1 / 3 . This scaling relation is reminiscent 
to that appearing in KP or KdV equations [35] . 

Of course this result is conjectural because the argument is based on hypothesis that are 
not fully supported by calculations: the use of the linear approximation near the shock, and the 
form of the dispersion relation. Nevertheless we find the same estimate using a quite different 
approach in the following section. 


3 Quantum Virasoro flows and quantum hydrodynamics 

In this short Section we would like to make a step towards connecting the previous hydrody¬ 
namic formalism to possible quantum flows on Virasoro module^]. The first step consists in 
reformulating the TT perturbations as random diffeomorphisms. This connection, based on the 
fact that the stress tensor is the generator of diffeomorphisms, is made here only to first order 
in g. We will use it to recover some of the previous results from the quantum flow perspective. 

3 This tentative formulation was motivated by a seminar given by J. Cardy in December 2014, at the conference 
“Mathematical physics of non-equilibrium quantum systems”, King’s College London, in which he presented a 
possible connection between TT perturbations and CFTs in random metrics [2Bj. The approach we follow is 
different from that of J. Cardy, as we use a purely algebraic setup without appealing to Hubbard-Stratonovich 
transformations. 
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To be more specific, let <fi(x) be a chiral primary field on the line - with scaling dimension h - 
and consider its time evolution <j)(x, t ) with respect to the perturbed hamiltonian (HD. Because T 
commutes with cj) and since T is the generator of diffeomorphisms, we claim that the hamiltonian 
flow on chiral fields such as cj) is that induced by the “vector field” v(x) := 1 + gT(x ), that is: 

</>(*, t) := e~ itH e +itH = (X[(x)) h c/>(X t (x)), (48) 


with X t the flow line generated by v and started at x at time t = 0, i.e. 

X t (x) = v(X t (x)) = 1 + gT(X t (x)), with X t=0 (x) = x. (49) 


The statement (|48j) is actually valid to first order in g only, because claiming that <p(x,t) com¬ 
mutes with T is true only to first order in g. 

To prove ea. (H8D we have to prove that (j>(x,t) = [X' t (x)] h 4>(X t (x)) is solution of dt<j)(x,t) = 
—i[H,(j)(x,t)\. Let us first recall well-known relations between diffeomorphisms, flow lines and 
the stress tensor. Since T is the generator of diffeomorphisms, its commutation relation with a 
primary field 4> of conformal dimension h read 

-*[ f dy£(y)T(y ), </>(x)] = h£'(x) <j>{x) + £(s) <t>'{x). 

Hence, for (f>(x,t) = [X' t ] h 4>(X t ) we have 


—i 


dyay)T(y),[X't] h ttX t ) 


hi\X t ) [X’] h <l>(X t ) + £(X t ) [X'] h <t>’(X t ), 


provided that Xt commutes with T. The latter commutator coincides with the time derivative 
dt<j>(x,t) provided we choose X t such that X t = £(Xt), because then X[ = X' t ^'{X t ) by the 
chain rule. That is: to solve for the time evolution we have to choose Xt as the flow line 
started at x associated to the vector held £. Now, our perturbed CFT hamiltonian is H = 
f dy{T{y) + T(y) + gT(y)T(y) s ). Since T commutes with any chiral held, we can apply the above 
remark to first order in g provided we identify the vector held £ with 1 + gT. 

Thus, at least to first order in g, the dynamics of chiral helds is encoded into properties 
of the hows ( 1491 ) that we would like to call quantum Virasoro hows. These hows are random 
because T is a huctuating data. We do not yet know whether these quantum hows can be made 
mathematically rigorous — of course defining the probability distribution functions of the how 
trajectories from data on the Virasoro algebra requires appropriate regularization. 

Let us nevertheless use this how to deduce properties that were found in the previous sections. 

The temperature dependent deformation of the light cone is easy to understand using quan¬ 
tum Virasoro hows. The light cone velocity is simply the mean velocity (1 + gT(Xt)), evaluated 
in the initial state. It is equal to iy r = 1 + 7) 2 r , with T^ r the left/right temperatures, 

depending in which region the shock front is moving. 

The quantum Virasoro hows also give an alternative understanding (based on scaling argu¬ 
ments) for the f 1//3 -scaling behavior we have been arguing for in the previous Section. Indeed 
let us consider two (quantum) trajectories started at neighborhood points x and x + 6x, and 
let 5Xt be their separation after a time duration t. By construction, its time evolution is 
(5Xt) = gST(Xt). Since T is of scaling dimension 2, 5T typically scales like l/(5Xt) 2 , or 
equivalently (6Xt) ~ g/(5Xt) 2 . By integration this yields that (6Xt) 3 scales linearly in time, 
(SXt) 3 ~ gt, which indicates that the correlation length in quantum Virasoro hows scales with 
time as (gt) 1 / 3 , which thus gives the natural width of the shocks. It would be interesting to 
make such arguments less conjectural. 
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4 Conclusion 


In this paper, we reviewed the general hydrodynamic approach to non-equilibrium quantum 
field theory. In order to study non-equilibrium effects it was important to keep two conserved 
quantities, including a parity-odd conserved charge, in the local fluid description. We then 
applied these ideas to develop a hydrodynamic approach for conformal field theory perturbed 
by the TT operator. We verified by direct quantum computations that the main aspects of 
hydrodynamics emerge at first nontrivial order in the coupling g. 

An immediate question is about the higher orders in perturbation theory. Clearly there were 
many simplifications at first order, in particular the emergence of Lorentz invariance, which are 
not present at higher orders. At higher orders, the fact that the perturbation is irrelevant may 
become more crucial and generate non-universal contributions. 

It will also be interesting to study the internal structure of the shocks, already at first order 
in g. We believe that the conjectured scaling relations are correct and universal. An important 
further question concerns the existence or not of universal structures inside the shocks, in a way 
similar to what has been revealed in classical fluctuating hydrodynamics [37]. 

In this paper we have focussed on mean transport phenomena and it will be very interesting 
to describe fluctuations and their large deviation functions. In particular one may wonder 
whether the extended fluctuation relations [38] . proved in pure non-equilibrium conformal field 
theory and in integrable systems (and conjectured to hold as well in higher-dimensional CFT 
m), remain valid once the TT perturbation has been turned on. 
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A Relativistic thermodynamics 

Consider a one-dimensional quantum system with the following general structure: a hamiltonian 
H generating time translations and a momentum operator P generating space translations, both 
expressed as integrals of local densities: 


H = 


dx h(x), 


P = 


d xp(x)\ 


and the conditions that (1) both are conserved: 


d t h + d x j = 0, d t p + d x k = 0, 


(50) 


(51) 


and 

tive 


that (2) the momentum density is equal to the energy current, plus possibly an even deriva- 
of the energy current: 


P = j + h dl e j, i > 1 


(52) 
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where /r is some constant. We are interested in the “generalized thermal” state, where both 
the hamiltonian and the momentum operator are present. We parametrize it by a “rest-frame 
temperature” T and a “boost parameter” 9: 

Tr ( cos h ^ — sinh 0 J 3 ) _ _ ^ 

^ ^ r X 1 r ^e — t( cos ^ 9 H— smhflp)^ 

and below we will denote h^ = (h)g, etc. 

We assume that there is parity symmetry whereby H is invariant and P changes sign, that 
the ground state is parity invariant, and that the limit T —>• 0 of the state (■ ■ -)g is the ground 
state, independently of 9. Further, we assume that the densities are homogeneous and (as 
the notation suggests) not explicitly time dependent: for instance, [P , h{x)\ = id x h(x ) and 
[H, h(x)] = —idth(x). Note that there is a simple gauge symmetry: we may change h(x) H > 
h(x) + a 1 and k(x) e-)- k(x) + a'l for any constants a and a' without changing any of the 
conditions. We will also use the principle according to which 

[P, O] = 0 => O oc 1 (53) 

for local densities O. Note that pg = jg since the average of a derivative of a homogeneous 
density is zero. 

It turns out that this structure is constraining enough to imply relativistic invariance. We will 
indeed show from the above that averages in a generalized thermal state take their relativistic 
form, and that a standard thermodynamic relation between the variations of pressure and of 
temperature hold: under the choice of gauge such that ho + ko —>• 0 as the temperature T goes 
to 0, 


hg = cosh 2 * 9 ho + sinh 2 9 ko 
k g = sinh 2 9 ho + cosh 2 9 ko 

p g = sinh 9 cosh 9 (ho + ko) (54) 


and 


(h 0 + k 0 )dT = Tdk 0 . 


(55) 


Remark that, from (l55lh if ho ~ AT q for some q > 1 as T —>• 0, then ko ~ A'T + AT g /(q — 1), 
where A and A' are some constant^] (the case q < 1 would pose physical problems as it implies 
that the pressure becomes negative in some temperature range). Equations (1541) can also be 
written in the Lorentz-invariant form T = ko : rf ;V + (ho + ko) UgUg, where r= diag( —1,1) is 

k ' 9 | is the stress-energy tensor and Un = ( cos | 1 ^ ) i s the 

Pe k g) 9 \smh9j 


the Minkowsky metric, V f) = 


relativistic velocity. 

These are purely dynamical relations. All the information about the particular model at 
hand may be imbedded into the “equation of state”, and once this information is given, the 
above provide all averages explicitly in terms of the temperature and the boost parameter. By 


4 In fact, there is a stronger statement, based on a weaker assumption: if ho < AT q V T > 0 small enough, 

then ko = A'T + g(T) with g(T) < -^A[T q VT>0 small enough. 
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parity invariance in the thermal state, po = 0. Hence there is one parameter characterizing ho 
and kg, so we can define the function F such that the following equation of state holds: 


k 0 = F( h 0 ). 


(56) 


Then, given F the relation (1551) provides the explicit temperature dependence in the thermal 
state, 


f*° <M _ f h0 d £F'{1) 

J i + F- 1 (£) ~ J £ + F(ey 


(57) 


Note that in the generalized thermal state, the equation of state is more complicated, and, from 
relations m, takes the form 


\J (he + kg) 2 — 4 P g 


he + k 0 = 2 F 



(h e + ke) 2 


4pe + he 



(58) 


Although we restricted the analysis to one-dimensional systems, the above hold as well 
in higher dimensions, where the effective one-dimensional system is obtained by averaging on 
hyper-surfaces transverse to the x direction in a d-dimensional space. Denote by the cur¬ 
rent associated with invariance under translation in the direction x u , with x° = f, x l = x and 
transverse direction x = (x 2 ,... ,x d ). Then h(x), p(x) and k(x) are transverse averages of the 
energy density H _1 f v d 1 xT llll (x), of the x -momentum density V^ 1 f v d^x T 01 (x), and of the 
x-pressure density H _1 f v d^x T n (x), respectively, where V is the transverse volume. This is 
particularly relevant for energy transport, as the direction of transport makes the system effec¬ 
tively one-dimensional. We note that, for instance, in scale invariant systems in d dimensions, 
the equation of state is 

k 0 = dh 0 (59) 

Proofs 

First we may redefine the energy density and current as follows: 

h = h + n dfh, j =j + n dfj. (60) 

This is such that h is still an energy density, H = f dx h up to local densities at infinity, that 
the conservation equation holds, dth + d x j = 0, and that the new energy current is exactly equal 
to the momentum density, j = p. Further, these re-definitions do not affect the properties of 
the energy density and current under a parity transformation. Hence without loss of generality, 
below we assume that h and j have been chosen in such a way that j = p. 

For lightness of notation we omit the index 9, except when it is equal to 0 (thermal state). 
First we show that there exists a boost operator B giving rise to the algebra 


[B,H] = P, [B,P] = H, [H,P] = 0. 

[B,H] = —^i J d xx[h(x),H] = j dxxdth(x) = — J d xxd x p{x) = J d xp(x) 


Indeed, consider 


Then we have 


(61) 

( 62 ) 


= P 
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and 


[B,P\ = —i j d x x[h(x), P] = — J d xxd x h{x) 


dx h(x) = H , 


up to local densities at infinity. Local densities at infinity do not contribute whenever the 
operators are exponentiated either when acting on local observables by adjoint action, or in 
some density matrix where averages of local observables are takerH. □ 

Remark that we used the fact that the energy current is the momentum density only in 
the first calculation: for any translation invariant system we have [B, P] = H , even without 
relativistic invariance. 

Naturally, the algebra (1611) implies that the generalized thermal state is a boosted state, so 
that 

O e = (e eB Oe- eB ) 0 


for any local observable O. 

We next show that 

d e h = 2p, d e p = h + k, dgk = 2p (63) 


under the choice of gauge characterized by the fact that ho + ko —>• 0 as T —>• 0. 

Let B(x) = e~ lPx B e lPx be the translated of B. We have B(x) = B + ixH. Let now 
[B,h](x) = e~ lPx [B,h]e lPx be the translated of [B,h] with h = h{ 0). Since [B,h](x) = 
[B(x),h(x)\ we have: 

[B,h](x) = [B,h(x)\ -xd x p(x), 

where we use the equation of motion i[H, h(x)\ = —d x p{x). Integrating, 


P=[B,H] 


= / dx (-p(x) + [ B , h](x)) 


-P + 


dx [B , h\(x) 


which is valid up to local densities at infinity, whereby 

[B, h](x) = 2 p(x) + id x b{x) 


for some local density b(x). 

Let us show that b(x) is a homogeneous density, d x b(x ) = —i[P, 6(x)]. By definition [B, h](x) 
is a homogeneous density: —i[P, [B, h](x)] = d x [B, h\(x). Since alsop(x) is, then d x b{x) is a local 
homogeneous density: d x [P,b(x)\ = id x b(x). The last equation means that [P, b(x)] — id x b(x) 
is independent of x. According to (1531) . the only local densities that are independent of x are 
those proportional to the identity operator 1, hence 

[. P , 6(x)] — id x b(x ) = a 1 (64) 

for some constant a. Further, under parity we have that [B,h](x) — 2 p(x) changes sign (in 
addition to x changing sign). Hence b(x) changes sign, again up to a local density that is 
independent of x. That is, under parity, b(x) i->- —b(—x) + dl for some other constant a. 
Putting this into (|64D . we find that the left-hand side changes sign under parity, whereby we 
must have a = 0. Therefore 

[B,h](x) = 2 p(x) + [P,b(x)]. (65) 

s That is, this is valid from the point of view of the properties of time and space translations or of states, not 
the point of view of the averages of H and P in a state. 
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and this leads to the first of (1631) . 

By a similar argument, 

[B,p(x)\ = xd x k(x) + [B,p\(x) 

and integrating, 

H = [B , P} = J dx (-fc(x) + [B,p](x)) 

up to local densities at infinity, therefore 

[B,p\(x) = h(x) + k(x) + id x b'(x ) (66) 

for some local density b'(x). Again d x b'(x ) is a homogeneous density, so that, by arguments as 
above, [P, b'{x )] — id x b'{x ) = a'\ for some constant a '. We cannot use parity symmetry in order 
to fix the constant a ', however we may absorb it into an appropriate gauge choice for the sum 
of the energy density h(x) and the pressure k(x). With this choice of gauge, the second of (fTHTI) 
follows. Since as T —> 0 the state specializes to the ground state, which is assumed to be Lorentz 
invariant, then from (1631) this choice of gauge can be characterized by the fact that h + k —>• 0 
as T —>• 0. 

Finally, let us set x = 0 in (|66l) . giving 

[B,p\ = h + k + [P, b'] — a 1. 

We apply [ H , •] on both sides and use the algebra (Ifill) as well as the conservation equations Q. 
We obtain on the left-hand side 

~[P,P] + [B, [■ H,p ]] = ~[P,p] + lB , [P, k ]] = [P, [B , k]-p] + [H, k} 
and on the right-hand side 

[H, h] + [H , k] + [P, [H, b ']] = [P,p] + [H , k] + [P, [. H , b% 

Comparing, and using (l53l) . we get the relation 

[B,k\=2p+[H,b'] + a"l 

for some constant a This gives dgk = 2p + a". Again, as T —>• 0 the state specializes to the 
ground state which is Lorentz invariance, whence 2po + a" = 0. Further, by parity invariance 

of the ground state po = 0, whereby a" = 0. This shows the third relation of (1631) . □ 

From (1631) . relations (1541) immediately follow: keeping the temperature T fixed, the system 
of three first-order ordinary differential equations (j63[) in 0 has a unique solution with the initial 
conditions ho, po and ko at 6 = 0, and one can check that (|54l) is a solution. □ 

Finally we show (1551) . For this purpose, we show the two relations 

sinh 9 dgp — cosh 8 dgh = sinh 9 T — cosh 9 Tdrp 

sinh 9 dgk — cosh 9 dgp = sinh 9 T c^p — cosh 9 T (9^k (67) 

Denote the density matrix as e~^ hH ~^ pP . These two relations are a re-writing of 

d d d d 

Wh V = W P h: = W P V 
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respectively. The first one is trivial as the equality boils down to f dx (h(x)p(0)) c = f dx (h(0)p(x)) c 
where (• • • ) c is the connected average. The second one in (1671) can be obtained by differentiating 
the first one with respect to 9 (assuming that the T and 9 derivatives can be interchanged), and 
by using (1U51) . Then consider the second relation of (HTTP at 9 = 0, which is <9e(p)|6»=o = Tdr{k) 0 - 
Using (1631) this is (h) 0 + (k) 0 = Tdr{ k) 0 . □ 
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